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ABSTRACT 


The effects of turbulence on the performance of imagers or 
on beam forming optical systems are well expressed by the 
optical transfer function (OTF) or its magnitude, the 
modulation transfer function (MTF). It has been shown that the 
MTF can be expressed in terms of the Fried model by means of a 
single number, the turbulence structure constants for optical 


index, C”, provided that a prorerly path weighted value is 


nN? 
obtained. Based on current theory the weighting needed in the 
path-position weighted value of on is a function of (z2/z,.)5/, 

heavily emphasizing the region near the imaging system or beam 


forming optics. It is the conciusion of this report, that the 


path weighting function (2/Z,)>/” is correct. 
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iil NTRODUGIRLON 


The quantity Oe which varies with the micrometeorology, 
can be determined in its properly path-weighted form with the 
slit scanning telescope system developed at NPS by Crittenden 
et al (Crittenden, 1978). ok can also be computed from 
measurements of the temperature structure constant, Cie Poanit 
measurements of of are based on assumptions of the statistical 
form of the turbulence, while optical determination of? C4 with 
the slit-scanner is independent of this statistical nature. 

Anomolies in the distribution of Ch due to air flow around 
obstructions near the optics is a common problem. fTh:.s has a 
direct naval application since for ship mounted optics, the 
heavily weighted part of the path is the shipboard end. 

The objective of this work is to design and conduct an 
experiment to determine the path-weighted values of Ck and 
compare these results to the existing model for the path 
weighting of the MTF. The model used is that proposed by the 
theoretical treatment of Fried (Fried, 1966) which predicts a 


(ey...) °/ Waependence where Zo = range. 





II. THEORETICAL CONSIDERATIONS 


A. GENERAL 

Turbulence and temperature gradients in the atmosphere 
result in random fluctuations in temperature about an average 
value. Similarly, turbulence and humidity gradients result in 
humidity fluctuations. The optical index of refraction depends 
On both the temperature and the humidity, so that the index, n, 
undergoes fluctuations around its mean value. 

A quantity useful in characterizing optical turbulence 
effects is the average of the square of the difference of the 
index of refraction at two points r andr separated by a 


distance r, i.-e.,: 


For isotropic turbulence, the quantity See expressed by 


Tatarski (Tatarski, 1961): 


ote ot gee et Lo 


<[n(r,) - Gies)| => = 2h Ci, £6 


is called the “index of refraction turbulence structure 
constant". This relationship holds within the "inertial 
sub-range", 

ee << r << Lo: 
where f and L, are the inner and outer scales of turbulence 
respectively. These scales correspond roughly to the size of 


the smallest eddies of turbulence, about a millimeter, and the 





largest eddies, often taken to be the order of the heignt above 
the surface of the point or propagation path of interest. 

The quantity C,° is a property at a point in the atmosphere. 
The mean of the integrated value of Co over the optical path, 
properly path weighted becomes the primary parameter of 
interest for the optical properties of the atmosphere due to 
turbulence. 

Values of Geeat specific points in the atmosphere can be 
measured from the temperature and humidity fluctuations. The 
index of refraction n, is related to the absolute temperature, 


ee DY 
m- 21=79 x 107 p/t 


with the pressure p in millibars. The average of the square of 
the temperature fluctuation can be shown to follow an equation 
Similar to that for index fluctuat:ion, leading to a definition 
of the temperature turbulence structure constant Cr in terms 


OF: 


cae) Sere > = Cie r?/? 
2 


Using the expression for index in terms of p and T above, this 


i 


leads to (Ochs, et al, 1969): 


Goel 79 me 1OT® (p/T> > ce 





Earlier work at NPS (Crittenden, 1978) indicates that the 
contribution to the index variation from humidity fluctuations 
is small compared to temperature fluctuations and is therefore 


neglected in this report. 


B. TURBULENCE EFFECTS ON OTF AND WANDER 

Turbulence in the atmosphere, as expressed by a 
path-integrated oe determines the performance of imaging or 
beam-forming optical systems, and this performance can be used 
to measure the path-integrated Cy. 

The effects of atmospheric turbulence on the performance of 
an imaging or beam projection system are usually treated 
(Fried, 1966; Lutomirski and Yura, 1974) by considering the 
atmosphere-plus-optics combination as a linear system, with the 
image being the response of the system to the object. If there 
were no sources of image degradation between the object and the 
image, a point object would have a point as its image. 

However, diffraction, caused by the finite aperture of the 
system, is always present. Within the present context, 
refraction and diffraction caused by the index of refraction 
inhomogeneities in the atmosphere are also present. These 
effects cause the ideal point response to a point object to be 
smeared out. If the point spread function has the same shape 
regardless of its position in the image plane, then the image 
function becomes the convolution of the object function and the 
point spread function. In this case the convolution theorem of 


Fourier transform theory yields: 


10 





LV) = HW yrvy) © OM xr Vy) 


where sai il H (wv amano and ait a are the two-dimensional 
Fourier transforms of the image function, point spread 
function, and object function respectively and the values are 
the spatial frequencies. The quantity H(vy,v ) is known as the 
optical transfer function (OTF) and its magnitude is called the 
modulation transfer function (MTF). The optical transfer 
function is thus the fundamental quantity which contains the 
effects of the atmopsheric turbulence. The modulation transfer 
function is also used to characterize the effects of 
turbulence, but it lacks information on phase. 

For isotropic turbulence and circularly symmetric optics, 
the OTF is a function of v,=vy=|v|. In this case, a 
one-dimensional treatment is adequate. Physically, a 
one-dimensional image signal is cbtained by scanning the image 
of a point object with a slit. Then the point spread function 
ls replaced by the line-spread function: 

h(x) = [ace dy and i(v) = H(v) oy) 

Theoretical studies (Fried, 1966; Lutomirski and Yura, 
1971) for an imaging system wnich include the effects of 
atmospheric turbulence show that the total optical transfer 
function can be expressed as the product of the transfer 
function for the optics and the transfer function for the 


atmosphere, 
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Heot (Y) = Hope (Y) * Hap (VY) 


Because of the stochastic nature of turbulence, quantities 
which relate to it can be understood only by means of an 
ensemble average. The problem of calculating the atmospheric 
MTF is considered in two limiting situations: the short term 
or short exposure (S.E.); and the long term or long exposure 
(L.E.). The short term MTF describes images obtained in the 
limit of time intervals sufficiently short that the turbulence 
can be considered frozen. The average short term MTF is then 
the average of the MTF's corresponding to many such images. 
The long term MTF results from taking an image of sufficiently 
long term which sees effectively all possible turbulence 
configurations. The average kong term MTF is, therefore, the 
Same as that from a single long time integrated exposure, since 
every turbulence configuration has been taken into account 
(properly weighted). 

If one observes a series of line-spread images, one would 
find that each image is broadened over the diffraction limit 
imposed by the optics and that the center of the area of each 
image appears to wander. The first effect is called image 
spread and the second image wander. The short exposure MTF is 
obtained by taking the average (and then the Fourier transform) 
of successive line-spread functions in which each line-spread 


function is shifted so that the center of the area of each 


a2 








curve has a conmon origin. This procedure removes the image 
wander. The long exposure MTF is obtained by not performing 
this shift in origin, therefore including image wander. The 
long exposure MTF is seen to be the result of image spread (the 
short term MTF) plus image wander. 

Several theoret.ical models have been proposed from which 
the L.E. MTF can be: derived. The only theoretical treatment of 
the S.E. MTF thus far is that of Fried (Fried, 1966). The 


result can be summerized by the following expression: 


M(f) = <H(f)> = exp | -57.64 B GA Zs y72/3 8/3 7 - (21h) | 
where 
XK = wavelength 
f = angular spatial frequency (cycles/radian), = fv 
Vv = linear spatial frequency (cycles/meter) 
F = focal itengt:h of system 
D = dianeter of! optics 
a = Lane 
f tor the L.E. OTF 
a =41/2 for the far-field S.E. MTF, D << (za)?/2 
1 for the near-field S.E. MTF, D > ee 
, 1 for a plane wave 


3/8 for a spherical wave 


Use of this quantity, f, the angular spatial frequency 


appearing in this expression, avoids the necessity to know the 


a 


Ol il —_— Ss 


— > Gantt bans 
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focal length or magnification of a system. Most experimental 
Situations which employ a diverging laser beam as the light 
source or object can be considered to correspond tc the 
spherical wave model. Plane waves are rarely enccuntered and 


will not be considered in this report. 


C. PATH POSITION WEIGHTING OF Cy, FOR THE MTF 

The expression (in the preceding section) for the MTF 
applies to the case where a 1s independent of the position 
along the propagation path. If Cy, varies along the path, this 


equation generalizes for the L.E. spherical-wave ‘iTF to 


Z 
Oo 
M(£) = exp [ -57.64 ee e5/> f ee) (2/29)? /? az | 
Oo 


That is, Cy, is replaced by a weighted mean cf 6given by 


The experimentally determined MTF can be used to infer this 
weighted mean for are 

The above equation expresses the path weighting of om 
indicated in Figure I-l in which the path region nearest the 
imaging, or beam projecting system at Z = Zo has the largest 
contribution. The opposite end of the path has zero 


contribution. 
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IiI. GENERAL PLAN OF EXPERIMENTAL PROGRAM 


A. GENERAL CONSIDERATIONS 

The general cbjective is to confirm the path-weighting 
model for the MTF. If the background es in a controlled 
environment is several orders of magnitude below the ee GE a 
highly localized temperature disturbance, then the optically 
measured ce will be due primarily to the injected disturbance 
and the contribution of the remainder of the path may be 
ignored. 

The major requirements of the experiment are then to 
construct an envi.ronment with sufficiently low and uniform 
eeeagrouna CC. and to construct a heat source which adequately 
meets the isotropic turbulence requirements. This heat source 
will require a method of monitoring to ensure a ake uniform 
disturbance as it 1s repositioned along the path in the 
controlled environment. If slight non-uniformities exist, they 


may be used as normalizing factors for the optically determined 


values of oe 


Bae PROCEDURE 

A sixty-one meter tunnel was constructed for the controlled 
environment and a suitable heat source designed for the 
temperature disturbance. By placing the heat source at 
progressive locations in the "tunnel", measuring the optically 
determined values of C, from the MTF as determined by the 
slit-scanner and on-line data processing system, the results 


16 





were plotted as oom versus the temperature disturbance location 
in the tunnel. These are compared to the theoretical curves 
fitted to the data taken on the same graphs. 


Monitoring of the heat source was achieved by using 


2 


platinum temperature probes across an AC bridge to measure Cy 


and thereby compute ore at the location of the temperature 
disturbance location. These values will not necessarily 
coincide with the optically determined values of Cf due to 
their dependence on the statistical nature of the temperature 
of fluctuations of the heat source (micrometeorology). The 
optical measurements assume a uniform or over the path. As 
mentioned earlier, the controlled environment ambient C4 is 


much less than the om of the heat source. 


Jeg 





IV. MEASUREMENT TECHNIQUES 


1. Long Term, Short Term, Wander 

The performance of an imaging system, or of a beam 
projection system such as the HEL and pointer (tracker) can be 
characterized and predicted in terms of the MTF of the 
atmosphere. The experimental program has been designed to 
provide verification of the theoretical model for the path- 
weighting of the MTF. 

As indicated in Section II, the MTF is expected to depend 
on the level of turbulence as expressed by the path-weighted 
value of So and on the wavelength and range. Two experiments 
were carried out to determine the effect of a temperature 
disturbance on the path-weighting of the MTF. These have been 
carried out in a controlled environment as free from other 
turbulence as practicable with uniform values of ee over the 
path except at the points of interest. 

Measurements of the MTF yield long term, short term, and 
wander. The internal consistency of these measurements are to 
pe compared. 

To calibrate the MTF in angular spatial frequency with a 
maximum scale of 500 microradians consisting of 1024 data 
points, a calibration signal is injected by placing a 
diffraction grating in front of the optics. The spacing 
between first order peaks is known to be 142.6 microradians, or 


A/d = 142.6/2 = 71.3 microradians. The dwell time is the 
18 





number of microseconds per point. For a dwell time of 20, the 
time between the first order interference peaks 


Zone eetisec/point) (1024 peints) (142.6/500) 


5.841 milliseconds 

Using a Hewlett-Packard 1743 oscilloscope, this time is set on 
the calibration signal by adjusting the gain of the scanning 
mirror voltage. Center display of the signal is obtained by 
adjusting the scanning mirror D.C. offset. 

The calibration grating is then removed and the desired 
signal observed. 

A Wavetek model 145 function generator provides the sweep 
Sagnad to the scanning mirror controller and a TTL output to an 
IEC P-12 pulse generator which provides a delayed trigger to 
the data processor NIC-80 described in section IV-2. This 
delayed trigger prevents the NIC-80 from seeing the sweep-back 
Signal from the scanning mirror. Figure IV-l diagrams the 
optical data processing equipment. 

2. On-Line Data Reduction 

The on-line data reduction unit is a Nicolet Instrument 
Company NIC-80 data processor. This is a data averager 
interfaced directly to Hewlett-Packard 9825B programmable 
calculator. The HP-9825B acts as the controller and 
input-output device to the NIC-80 and controls the plotter and 
printer. 
The NIC-80 digitizes the analog signal from the 


detector, processes the raw data to obtain the experimental 
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MTF, wander and ope The processed experimental data is then 
transferred to the HP-9825B which is better suited for floating 
point calculations. The HP~-9825B calculates the best fit curve 
and plots out the results in final form. 

3s Optical Equipment 

The optical equipment consists of an 8 inch diameter 
Celestron, Schmidt telescope with an equivalent focal length of 
160 inches and an overall length of about 50 centimeters. It 
has a .344 central obscuration ratio. Mated to the image end 
of the telescope is the scanning mirror assembly. It consits 
of an optically flat mirror mounted on a servo motor to provide 
the sawtooth scan of the image across an adjustable slit. The 
mirror turns the optical path 90 degrees. The scan produces 
the line spread function as a function of time. A trigger 
signal associated with initiation of each scan indexes the 
position of the image. 

The detector, mounted after the slit in the 
slit-scanner, is an RCA silicon photodiode type C30872. It is 
Operatec at 370 volts D.C. with a 9 volt battery powering the 
source follower. The detector used does not require cyroqdenic 
cooling which is a large improvement in reliability and ease of 
use Over a detector that did during earlier trials. 

The telescope optics are supported on a rigid tripod 
base, weighted with 50 pounds of lead which prevents mechanical 
vibration of the telescope system. The optics and detector are 


sealed as a part of the tunnel during optical measurements. 
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The signal from the detector is preamplified by a 
high-gain, low-noise PAR 113 amplifier. The signal is 
processed immediately with the on-line equipment. 

The source for the MTF system is a Helium-Neon laser 
transmitted from the far end of the tunnel. The laser beam 
passes through two polarizers to reduce its intensity and is 
diverged by a 10 centimeter lens to provide a spot size at the 
MTF telescope of about .5 meters wide. This provides uniform 
illumination over the aperture of the telescope and closely 
approximates light diverging as a spherical wave fiom a point 
source. The MTF telescope then images this pcint source. The 
laser is operated in the CW mode for this purpose. The laser 
wavelength utilized for MTF measurements is the red HeNe at 
0.6328 micrometers. 

4. Controlled Environment (Tunnel) 

The tunnel was constructed of three-eighths inch 
construction plywood with completed sections measuring 2 feet 
by 2 feet by 8 feet (61 cm by 61 cm by 2.44 m). 25 tunnel 
sections were made for a tunnel length of 200 feet or 61 
meters. Prior to construction an underseal of water repellant 
was applied both inside and out. The inside was then painted 
with flat black and the outside with an exterior white. The 
individual sections were constructed with 12 metal edge 
brackets, 8 end corner brackets, and 56 wood screws. The 
joints were sealed with a heavy duty tape to prevent air flow 


through any of the seams. 
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A special heat source tunnel section was built by 
modifying an existing section. An opening about 1 1/2 times 
the area of the heat source was cut in the top at one end anda 
64 cm high chimney constructed over it. A wooden baffle sheet 
was placed over the top of the chimney (supported by small 
blocks) in a manner to allow air flow but prevent down drafts 
from entering the tunnel. A wide slit was cut in the bottom of 
the tunnel section under the chimney to allow for air flow 
across the heat source by convection. The slit was baffled to 
prevent entry of side drafts. Finally a small slit was cut in 
the side of the tunnel section for the insertion of the 
platinum temperature probes above the heat source. The heat 
source was then placed under the chimney in the bottom of the 
tunnel section (Figure IV-2). 

The tunnel, located in the basement of Spanagel Hall at 
NPS, was assembled for the experimental runs only after working 
hours since it blocks several doorways and was disassembled 
after each run. When assembled, the joinina and sealing of 
tunnel sections was achieved by taping. During the experi- 
mental runs, the heat source was injected every other tunnel 
section or every 4.88 meters along the 61 meter path. 

A table was positioned at each end of the tunnel to 
house the optics and detector at one end and the laser and 
lenses at the other. During the runs, these tables were sealed 
and became an integral part of the tunnel's controlled 


environment. 
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If the tunnel were a perfect controlled environment 
with no temperature disturbances, a point source at one end 
would be imaged by the optics as a point source at the other 
end and the MTF would be a straight line (normalized Ch = 1 
independent of spatial frequency). Figure IV-3 shows the very 
high degree of "perfection" achieved for the controlled 
environment inside the tunnel. 

5. Heat Source 

As discussed in Section II, a region of isotropic 
turbulence for 2, <<r<< Lo is required (or approximately so) 
from whatever heat source would be used. In this region Cy 
must be reasonably uniform. For several heat sources the C as 
a function of r (the probe spacing for the or from Ce equipment. 
described in Sections V-B and C) was plotted. The heat section 
used gave the best results and the data is displaved in Figure 
Iv-4. 

The heat source consists of a two conductor heating 
element bent in reversing V's mounted on a wooden stand with 
ceramic standoffs. The wide spacings of the elements were to 
assist in providing the larger turbulence size which was 
considered to be about the height of the probes above the heat 
source, 22 cm (Lo). A screen of wire spacing .15 cm was placed 
over the neating element to provide the smaller turbulence size 
a). 

POWer was supplied by a Variac set at 70 volts AC which 
produced a temperature 22 centimeters above the heat source of 


about 27-30 Celsius. 
23 





Assuming £.<<Yo9<<Lo or .15 cm <<r<< 22 cm the region 
of isotropic turbulence would be expected to be in the vicinity 
of 1.5 to 2.5 centimeters. Figure IV-4 shows a peak around 2.0 
centimeters the remainder of the curve falling off as r'/3 (the 
path dependence of C,,). The 2 centimeter region was considered 
to be approximately isotropic. A probe spacing of 6 
centimeters was eventually decided upon due to the sharp peak 
at 2 centimeters. A more linear section of the curve is 
desired for more uniform results since the measured oF will be 
used to monitor the stability of the heat source and provide 
weighting for the optically determined results. 

Thus decided upon, the heat source was placed in the 
bottom of the special tunnel section designed for it (Figures 


IV - 5 and 6). 


2 
Be. Cy POINT DETERMINATION 


To determine ofp 


Car = <[T(r,) oe.) ae 


T must be measured. The value of r, the temperature probe 
spacing in meters has been picked at .06 meters as discussed in 
the previous section. 

Using two platinum wire temperature probes, ( %= .0036° K- _ 
each with a resistance of approximately 56 ohms, across an AC 


Bridge, the output is cailbrated for AT. The bridge 


sensitivity was set and measured at 6.67 volts/ohm. 
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1 
AT per volt (.0036 °K-") (56 OHMS) (6.67 Volt/OHM) 


AT per volt ~-/44 K per volt. 

Using a spectrum analyzer on the AT fluctuations, the half 
power point for the AT fluctuations was determined to be about 
& hertz. A Hewlett-Packard 3438 A System Voltmeter, sampling 
at an 1l hertz rate for 500 samples (45 seconds), is then used 
to determine the values of AT. The system voltmeter is 
controlled by a Hewlett-Packard 9825B programmable calculator. 
It squares each AT sample and averages the 500 squared samples. 


Multiplying by rh , the calculator outputs Ge. 


Z 
Soe, FROM Cc, 


2 
To determine Cy , 


amo 105° (P/T?) ) C2 


the absolute temperature of the heat source must be measured. 

Using a Copper-Constantan thermocouple referenced to 0” 
Celsius (ice water) and sampling at a 3 hertz rate using a 
Hewlett-Packard 3490A multimeter controlled by the 
Hewlett-Packard 9825B calculator, forty samples of temperature 
are taken and averaged. Twenty are before the sampling for ge 
and twenty after. 

T = <Thermo-Couple Voltage> 32809.524 ue per volt. Adding 


273.15 gives the temperature in degrees Kelvin. 
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Assuming a constant atmospheric pressure of 1013.25 
miliibars, the H/P 9825B calculates or to be used for 
monitorince the heat source and weighting of the optically 
determinec. values of CA. 

Figure IV-7 shows a block diagram of the equipment used 


in the Cf , T and Ca measurements and calculations. 
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V. RESULTS 


A. EXPERIMENTAL DATA 

m 22 March 1983 

During the evening of 22 March 1983, the first experimental 
run was made in the basement of Spanagel Hall at NPS. A 
“guiet" tunnel result similar to Figure IV-3 was obtained 
confirming the quality of the tunnel's environment. The point 
values of os from on were fairly uniform along the path with 
an apparent systematic increase on the detector end of the 
tunnel. The short term results (Figure V-l) were in excellent 
agreement with theory but the long term results (Figure V-2) 
though good, were not as convincing. The results from wander 
(omitted) were totally random in appearance. 

Since wander is the shift in image center, vibrations could 
cause these erratic results. The initial telescope mount was 
not the one described in Section IV-A, 2. The original 
mounting was noted as heing susceptible to vibration and the 
wander data on this run motivated its alteration. 

22) ees i 1983 

After the telescope vibration problem was corrected, a 
second experimental run was conducted during the evening of ll 
April 1983 in the same location. 

cet from e, point values (Figure V-3) were very uniform 
over the path. Again the short term results (Figure V-4) were 


in excellent agreement with theory while the long term results 
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(Figure V-5) were not as good, but still showed the general 
shape of the curve. The wander results (Figure V-6), though 
vastly improved, cannot be used to substantiate the theoretical 


model. 


B. ANALYSIS OF DATA 
The following conclusions are drawn from the experimental 
Gata : 
1. The short term results confirm the Fried model for the 
path-weighting of the MTF. 
2. The long term results give good credinility to the Fried 
model, but cannot be used to positively confirm it. 
3. The wander results do not support the theoretical model. 
The overall conclusions is that the Fried model is correct 
and that systematic noise is responsible for the poor quality 
of the wander results and aberrations in the long term results. 
As described earlier, wander is not removed during the 
processing of the long term data so any physical or electrical 
vibration of the experimental system will affect the wander 


primarily and the long term to a lesser degree. 


Cs. RECOMMENDATIONS 

Since lead bricks weighted all physical components of the 
system, the likelihood of mechanical vibrations during the 
second run is low. 

A large source of possibile error seems to be in the sweep 


generation of the scanning mirror. The flatness of the 
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sawtooth sweev voltage was adjusted by eye. To ensure 


linearity over the sampling region of the sweep, 
transform of the ramp function could be computed 
to a spectrum analysis of the scan-mirror output 
deviation from the predicted components would be 


mirror vibration and hence, wander. 
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